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Abstract 



In the past year, polarized deep inelastic scattering experiments at CERN and 
SLAC have obtained structure function measurements off proton, neutron and 
5^ \ deuteron targets at a level of precision never before achieved. The measurements 

can be used to test the Bjorken and Ellis-Jaffe sum rules, and also to obtain in- 
formation on the parton distributions in polarized nucleons. We perform a global 
leading-order QCD fit to the proton deep inelastic data in order to extract the spin- 
dependent parton distributions. By using parametric forms which are consistent 
with theoretical expectations at large and small x, we find that the quark distri- 
butions are now rather well constrained. We assume that there is no significant 
intrinsic polarization of the strange quark sea. The data are then consistent with a 
modest amount of the proton's spin carried by the gluon, although the shape of the 
gluon distribution is not well constrained, and several qualitatively different shapes 
are suggested. The spin-dependent distributions we obtain can be used as input 
to phenomenological studies for future polarized hadron-hadron and lepton-hadron 
colliders. 



1 Introduction 



In the last few years there has been a resurgence of interest in the spin structure of 
the nucleon. This was largely initiated by the measurements of the polarized structure 
function Qi by the SLAG- Yale [Q and EMC [@] collaborations. In the 'naive' parton model, 
Qi can, like the unpolarized structure function Fi, be expressed in terms of the probability 
distributions for finding quarks with spin parallel or antiparallel to the longitudinally 
polarized parent proton: 

Fi(x,g2) = |Een^(^) + ?»] (1) 

g 

9i{^,Q') = I E e^M^) + ^^1^)] ' (2) 

where 

? = 9t + ?i > M = (l^ -Ql ■ (3) 

The renewed interest in the subject was triggered by the first precision measurement of 
the integral of by the EMC collaboration 0, 

r?(EMC/SLAC) = /^rfa;^?(x,g2) =0.126 ±0.010 ±0.015 (Q^ = 10 GeV^) , (4) 
Jo 

which was significantly lower than the 'Ellis- Jaffe sum rule' value of 0.18 [Q. This latter 
prediction is obtained by assuming that the net contribution of strange quarks to the 
proton spin is negligible. In the context of this model, and rid, the net spin carried 
by up and down quarks respectively, can be obtained from the /?-decay rates of the octet 
hyperons: 

r]u = dxAu{x) = 2F 
Jo 

Vd= f dxAd{x) = F- D . (5) 
Jo 

Attributing the difference between the Ellis- Jaffe prediction and the SLAC-EMC mea- 
surement to a non-zero strange sea polarization leads to the rather unusual result that 
the net spin carried by the quarks, 

r]T.= I dx AS(a;) = f dx [Au + Ad + As] , (6) 
Jo Jo 

is very small. This so-called 'spin crisis' precipitated an enormous amount of theoretical 
discussion - a clear and comprehensive review can be found in Ref . . 
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One of the most compelling explanations for the violation of the Ellis- Jaffe sum rule 
is that a substantial amount of the proton's spin is carried by gluons. As first pointed 
out in 0, the polarized gluon contributes to gi via the 75-triangle anomaly. Thus in 
perturbative QCD, the naive parton model result of Eq. (H) is replaced by 

9i{^.Q') = \T.^I t- [^q{x/y.Q^)+Aq{x/y,Q' 
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+ \f/-^ AGix/y, Q') |n, ^ ACc(y) + . . ] (7) 



where AG = — G| is the polarized gluon distribution. Now according to the Altarelli- 
Parisi evolution equations 

VM') = + O{al) (8) 



rflogQ2 

a.{Q')VG{Q') = + O(a^), (9) 



lilogQ 

which implies that the gluon contribution to gi is formally of the same order as the 
quark contribution. The splitting of the structure function into two conserved quantities 
associated with quark and gluon contributions to the net spin is natural in the sense 
that it allows the former to be identified (up to finite mass corrections) with the SU(6) 
"constituent" quarks. Retaining only the leading terms in Eq. (|^), assuming three quark 
flavours, and using a factorization scheme where ACciy) = — y) (see Section 2), 

gives 

gli^.Q') = \ H e^fAg(x,Q^)+Ag(a:,Q2)l-^^^AG(x,g2), (10) 



2 

q=u,a,! 



- I dx g, - -vu + -v<i + Y^vs - -^-^nc (11) 

It is now straightforward to calculate the amount of net gluon spin needed to explain 
the SLAC-EMC sum-rule data. For example, with = and 7]^-, Tj^ again obtained from 
hyperon decays, we find r^c ~ 5 at ~ 10 GeV^ [0]. Although this might be considered 
a surprisingly large number there is no violation of spin conservation, since there is also 
a contribution to the proton's spin from the orbital angular momentum of the partons: 

\ = \ri^ + VG+{L.) . (12) 
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According to Eq. @, the second and third terms on the right-hand side increase as logQ^ 
in such a way that their sum is constant. 

In Ref. - hereafter referred to as AS - a simple leading-order QCD model, with no 
polarized strange sea at Q'^ = Ql = A GeV^, was used to extract polarized quark and gluon 
distributions from the SLAG and EMC data on gi{x, Q^). Although the integrated parton 
distributions were reasonably well determined, the shapes of the distributions, especially 
those of the ci-quark and the gluon, were poorly constrained. As a result, counting-rule 
and Regge arguments for the large and small x behaviour had to be invoked. The main 
purposes of the exercise were (i) to show that a consistent set of distributions could be 
derived, and (ii) to present parton distributions which could be used for polarized lepton- 
hadron and hadron-hadron collision phenomenology. 

Since the AS analysis was performed, there has been a dramatic increase in the amount 
of polarized structure function data available. Using a polarized ^He target, the E142 
collaboration at SLAC has measured the neutron structure function g" [||. The SMC 
collaboration at CERN first measured the deuteron structure function gf [P] and more 
recently has measured gf |10], improving on the earlier SLAC-EMC measurements. These 



new results allow a test of another important sum rule due to Bjorken |jTT| 



Bj 



dx ((?! (x, 



<(x,g2)) = ^ 

9v 



71 



+ 



(13) 



This sum rule follows from SU(2) isospin invariance and is a rigorous prediction of QCD. 
Any disagreement between the measured and predicted values would invalidate the QCD- 
improved parton model approach on which the present study is based. The most recent 
experimental measurement - using all available data JlOl - gives 



1^7 = 0.163 ±0.017 



= 5 GeV') 

to be compared to the third-order QCD prediction |]12 

(g2 = 5 GeV') 



r*j5 = 0.185 ±0.004 



(14) 



(15) 



The agreement is acceptable, especially considering that higher-twist contributions might 
still give a small contribution at this Q^. In our parton distribution analysis, therefore, 
we shall use isospin symmetry to relate the distributions in the proton to those in the 
neutron. 

It is interesting that the new measurement of the integrated proton structure function 
||10||, again using all available data. 



r?(SMC/EMC/SLAC) = 0.142 ± 0.008 ± 0.011 (Q^ = 10 GeV^ 



(16) 
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is larger than the old SLAC-EMC result. This means that our polarized gluon distribu- 
tion (which accounts for the difference between the measured value and the Ellis- Jaffe 
prediction) will now be somewhat smaller than previous estimates: in fact we shall show 
below that^ ?7g ~ 2, compared to the AS result r/c ~ 5. 

The goal of the present study is to use the latest available data to update and improve 
the AS parton distribution analysis. Following AS, we adopt the point of view that 
there should be no significant polarization of the strange quark sea. Although the sum- 
rule measurements give an indication of the size of 7]g, there is a large arbitrariness in 
the shape of the gluon distribution. We will explore several qualitatively different shapes 
allowed by the data. The outcome of the analysis will again be a set of distributions which 
will provide benchmarks for future deep inelastic scattering experiments and which will 
be useful for phenomenological analyses of other types of polarized scattering experiment. 

As we have already indicated, the first moments of our parton distributions are con- 
strained by hyperon decay and sum-rule data. In the remainder of this section, we describe 
the input data which we use, and the corresponding first moments which obtain. 

The structure function Qi occurs in the antisymmetric part of the hadronic tensor 
in deep inelastic scattering, which can be expressed as the matrix element of the axial 
vector current between two proton states |]1T|. These matrix elements can be further 



decomposed into the SU(3)/ singlet and octet pieces. In order to compare the sum rules 
with experiment and to perform a evolution of the polarized parton densities, we 
rewrite Eq. (0) as 

g^, = ^Aq,ix, Q') + ^Aqsix, Q') + JaS(x, Q') - lM^AG(x, Q') (17) 
iz oo y o Ztt 

with similar expressions for the neutron and deuteron: 

9i = -^Aqsix,Q') + ^Aqsix,Q') + ^AJ:{x,Q')-^^^AG{x,Q') (18) 

gt = ^Aqs{x,Q') + lAE{x,Q')-\^^AG{x,Q'). (19) 

Here we have defined octet and singlet quark distributions: 

Aqs = Au- Ad 

Agg = Au + Ad- 2As 

AS = Au + Ad + As. (20) 
The first moments of (p!7D -(|l9|) are: 

Ti = /3 + /s + /q - - Vg 



^ using also improved F and D values 
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(21) 



The axial- vector current matrix elements, J3 and Jg, are related to the F and D couplings: 

1 



12 
1 

36 



(F + D) 
(3F - D) 



(22) 



Assuming rjs = gives Jq = 4/8, leading to separate predictions for the quark contributions 
to the Ti's [|]. 

Measuring F and D in hyperon /3-decays gives their values at a momentum transfer 
scale of ~ (9(0.5 GeV^). However F and Z), being related to the octet axial vector 
current matrix elements, are dependent. To fix rju and r]d for our fits, we incorporate 
these perturbative corrections by replacing F and D by F and D, which are obtained by 
comparing the parton model expressions for and Jg + Jq with their values corrected to 
first order in perturbative QCD. Demanding that F and D reproduce the corrected values 
of I3 and Jg + Jo, one has, following IT] 



— (F + D) 
12^ ^ 



1 _ ^ 

TT 



36 



This leads to 



(3F - D) 

F : 
D : 



36 



(3F - D) 



(1 



TT 



+ 4(1- 



37r' 



(23) 



13q;s 
IStt 



ISvr 
2a. 



(24) 



Taking n/ = 3 and A = 177 MeV (see below), we calculate F and D at the starting 
point of our perturbative evolution Ql = A GeV^ (as(4 GeV^) = 0.2879), taking the 
experimental values for F and D from Ref. [p!3|: 



F) 



0.424 ±0.008 
0.718 ±0.007 



which fixes the first moments of the polarized u- and ci-quark distributions: 



Vu 

Vd 



0.848 ±0.016 
-0.294 ±0.011. 



(25) 



(26) 
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Finally, comparing the parton model expression for ([TT|) with the experimental value 
(0) gives 

r^G = 1.971 ±0.929 . (27) 
The above result for the moments and the sum rules are summarized in Table (|l|). 



partons 




0.848 ±0.016 


Vd 


-0.294 ±0.011 


Vs 


0.000 




0.554 ±0.019 


Vg 


1.971 ±0.929 


9i 

(4 GeV') 


r? 


0.142 ±0.015 


J- 1 


-0.048 ±0.015 


rf 


0.047 ±0.015 




0.190 ±0.002 



Table 1: First moments of the quark and gluon parton distributions and of gi 



The remainder of the paper is organized as follows. In the next section we discuss the 
different possible definitions of the polarized gluon distribution. In Section 3 we describe 
our fit to the polarized structure function data, in Section 4 we discuss the implications 
of our fits for the polarized sea quark distributions at higher Q^, and in Section 5 we 
summarize our results and present our conclusions. 



2 Definition of the polarized gluon distribution 



The form of ACg'(x) in Eq. (0) has been a matter of some dispute during the past few 
years [^, |T4|, |T^. Kodaira [|I^ first pointed out that the anomalous dimension of the 
singlet axial vector current was non-zero. According to Ref. |l^, the first moment of 
ACc.(x) is 



dx /S.Cg{x) 



-1. 



(2^ 



Based on these ideas, Altarelli and Ross ||^ evaluated the anomalous gluonic contribution 
io gi. A crucial point in this calculation is the regularization of coUinear singularities (from 
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g qq) in the photon-gluon fusion process. Depending on the regularization procedure 
used, one obtains different expressions for /S.Cg{x). In Ref. |^ the physical quark mass 
was used as a regulator, giving 

ACg(x) = (2x-l)lni^. (29) 

X 

The scheme dependence of ACg{x) was further discussed by Bodwin and Qiu They 
obtained the surprising result that, in dimensional regularization and for finite quark 
masses, the first moment of ACg{x) vanishes. Only the (unphysical) approach of assum- 
ing a small but finite gluon mass reproduced the result of [1^. A careful reanalysis [15| of 



the polarized photon-gluon fusion process showed that in fact this disagreement between 
different regularization schemes comes from the soft non-perturbative region. Contribu- 
tions from this region should be absorbed into the light quark density, and should therefore 
not contribute explicitly to /S.Cg{x). This leaves us with two different possible ACg(2;), 
both with first moment equal to —1. The first is obtained by introducing a non-zero 
gluon mass, while the second derives from either introducing finite quark masses or from 
applying dimensional regularization, always subtracting the soft contribution: 

ACg(x) (P < 0) = (2x - 1) {\n ^ - 2) , (30) 
ACg(x) (DR, mV 0) = (2x - 1) {\n - 1^ . (31) 

The gluon contribution to gi, Eq. (|^), involves the convolution of ACg{x) with AG. 
Rather than tying ourselves to a particular theoretical scheme, we adopt the AS procedure 
and define a "(71-scheme" gluon with ^Cg{x) = —6{l — x). Since only the first moments 
of the polarized quark and gluon densities are precisely defined, we are free to fix the 
higher moments in this way to obtain the simplest form for the gluon term in Eq. (|^. 
Of course, in a consistent next-to-leading order analysis, any change in the definition of 
ACg is compensated by a corresponding change in the 0{as) correction to the quark 
contribution, Eq. (|^), in such a way that the physical quantity gi is scheme independent. 
A different choice of ACg will then yield a different gluon and a different set of quark 
distributions. We are aware of the fact that our mixture of leading-order quark and next- 
to-leading order gluon contributions in gi is ambiguous. Nevertheless we consider the 
incorporation of the anomalous gluon, although formally of order a^, to be motivated by 
the scaling behaviour of a^AG and by the magnitude of this contribution to gi. When 
complete next-to-leading order corrections become available it will of course be necessary 
to perform calculations of different processes within a consistent scheme. 

Within our leading-order approach, it is straightforward to relate gluon distributions 
defined according to different schemes. Suppose that one uses a different coefficient func- 
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2 Gluon shapes 
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Figure 1: AG in different regularization schemes reproducing our (AS) set A gluon at Q^. 



tion, ACg{x) say, in the definition of gi. Then the corresponding polarized gluon distri- 
bution AG{x) (i.e. the distribution which reproduces the same gi at Qq) is related to 
ours by 

AGix, Ql) = - C ^AG (-, Qg) ACaiy). (32) 
This equation can be inverted in Mellin moment space. Applying a Mellin transformation 



to Eq. (|32D and defining 



AG{n) = dx x'^-^AG{x,Q 
Jo 

ACcin) = 



we obtain 



dx x^'-^ACcix) 



Zm Ja-ioo ACcyn) 



(33) 
(34) 

(35) 



For reference, we list in Table ^ the Mellin moments of the coefficient functions of Eqs. (|30| ) 
and (|31|). 

In the following section, we will determine a polarized gluon distribution at Ql from 
a fit to the structure function data. The gluon defined in any other scheme can then be 
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obtained by numerical integration of Eq. (^51). As an example, Figure |l| shows our gluon 
set A (see below) together with the two distributions which correspond to the different 
coefficient functions defined in Eqs. (^) and (pll). 







Mellin transformation ACcin) 


-6{l-x) 


-f 


-1 


(AS assumption) 






(2a;-l)(ln^-2) 
{e < 0) 


-1 


n'^ (n+l)'^ 


(2x-l)(lni^-l) 
(DR, ^ 0) 


-f 


n^+n2+n+l+(n^-n){7E+<I'(n+l)) 



Table 2: The different gluon coefficient functions and their Mellin transformations 



3 Fit to the polarized structure function data 
3.1 Parametric form of the fit 

As a result of over twenty years of experiments, the unpolarized parton distributions in 
the proton are very precisely known, see for example Ref. Ijl^. Some twenty independent 
parameters are needed to specify the various quark and gluon distributions at Q^. In 
comparison, measurements of the polarized distributions are still rather imprecise. It is 
therefore necessary to build in to the parametrizations a significant amount of theoretical 
prejudice. As measurements improve in the future, these expectations will be tested and 
refinements can be made. 

At present, the only constraints on the distributions are (i) the specification of the first 
moments of the distributions, as described in the Introduction and summarized in Table [^, 
and (ii) the requirements of positivity of the individual spin-parallel and spin-antiparallel 
distributions: 

/n = ^(/±A/)>0 ^|A/|</, (/ = g,G). (36) 

Therefore, the absolute value of a polarized parton disrtibution always has to be smaller 
than the corresponding unpolarized distribution. We therefore require a consistent set 
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of leading-order unpolarized distributions to provide the bounds - those of Ref. [|l^ are 
ideal for this purpose. The starting distributions are, at Q'^ = Ql = A GeV^, 

xd, = N,x'-'''\l-xf'''' 

xu = xd = xs = 0.1515(1 - x)'^-^'^^ 

xc = 

xG = 3.0170(1 -x)™ (37) 

and the evolution is performed with A = Alq = 177 MeV. The parameters N^d and A^^^ 
are fixed by the requirement that /q^ dx = 2 Jq dx d^ = 2. 



For consistency, we choose the same A and Qq values as |T8|, and similar starting 
parametrizations. Only our choice of performing the whole evolution with three quark 
fiavours differs from |T^. As explained in the Introduction, we start at Qo with valence- 



like forms for the up and down quarks and with an unpolarized sea: 

xAu^ = riuAux'^^il- xf^{l + jux) 
xAd, = r]dAdx''''{l-xp{l + -fdx) 
xAu = xAd = xAs = 
xAc = 

xAG = tigAgx"':' {1 - xf"^ {I + -fcx) (38) 
with normalization factors Aj {f = q,G) to ensure that Jq dx Af{x, Ql) = rjf: 

Non-zero polarized SU(3) symmetric sea quark distributions are generated dynamically 
for Q > Qo (see below). Note that (i) the positivity constraints (|36D demand that the 
bf parameter values are at least as big as the corresponding unpolarized parameters of 
Eq. (pTl), and (ii) the convergence of the first moment integrals requires af > 0. 

Even with these constraints, there are still too many free parameters in ( P^ ) for the 
amount of data available. We therefore impose additional theoretical constraints at large 
and small x, based on dimensional counting and coherence /Regge arguments respectively 
||19| , |20| , |2T| , p^ , . A good discussion of these can be found in Ref. . In the present 
context, they can be summarized as follows: 

(i) The assumption that the small x behaviour of the quark distributions is controlled 
by Regge behaviour requires = ad- 
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(ii) Coherence arguments at small x give oFq = a™^° + 1 for the gluon distribution. 

(iii) Spectator quark counting rules for the x ^ 1 behaviour suggest = + I and 
require b^ and b^ to be the same as in the unpolarized case. 

We have already mentioned that the most accurate data at present are for Qi. This means 
that the detailed shape of the d-quark and gluon distribution cannot yet be determined. 
Thus for the former we impose ju = 7d, and for the latter we suggest three qualitatively 
different forms, defined by 



Gluon A: 


7g = 


0.0 


Gluon B: 


1G = 


18.0 


Gluon C: 


IG = 


-3.5 



The reasoning behind these choices is as follows: 7g = (the 'standard' form) will give a 
fit with ~ Gi at large x, which corresponds to a fast-moving gluon being equally likely 
to have its spin parallel or antiparallel to the parent proton; 7^ large and positive will lead 
to G^ ^ Gi at large x, i.e. a fast-moving gluon which is preferentially polarized in the 
direction of the proton's spin; large and negative will give a large-x gluon whose spin 
is anticorrelated with that of the proton. For the latter two cases, we limit the absolute 
size of the parameter 7^ to ensure positivity of both G^ and G^, and in each case we fit 
for the remaining parameter bg. 

Table |^ lists the parameters which are fixed by the above theoretical considerations. 
In Section ^]3| we describe how the remaining parameters a^, •ju and be are determined 
from the data. 



Parameter 


Value 


Comments 


bu 


3.64 


unpolarized 6n|18l 


hd 


4.64 


unpolarized 6^1 18| 


aa 


1.000 


coherence arguments 21 


Vu 


0.848 


Vu = '2F 


Vd 


-0.294 


Vd = F-D 


Vg 


1.971 


from 



Table 3: Parameters fixed by theory 
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3.2 Evolution of the parton distributions 

Having fixed the starting distributions at Ql = A GeV^, the distributions at higher 
are obtained from the leading-order evolution equations M: 



«(Q2) ridy 
2n Jx y 



dAE(a;,Q2) 
(ilnQ2 

dAG{x,Q^) 
(ilnQ2 

dAqyai{x,Q'^) 

dlnQ'^ 2n Jx y 

where the splitting functions are 



AG{y,Q')AP,G 



X 



2tt Jx y 



X 



X 



AS(y, Q')APg, - + AG{y, Q')APgg - 



Aq,ai{y,Q^)AP^ 



X 



AP,G 
APGg 

APgg 



A l-{l-zf 
3 ~z 



{l-zf 



z (l + z) + ) z 
With the first moments AP^^ and AP^q both being equal zero it follows that 



(41) 



(42) 



dxAqj,ai{x,Q^' 



and 



dxAT.(x, 



(43) 



are both independent of Q^. With APq^ = 2 and APqq = Po/2, one finds, using the 
renormalization group equation, 



d 



dlnQ- 



-MQ^)'1g{Q^) = 0{al) where r/G(Q') = dxAG{x, 



(44) 



and so Ff, and F^ defined in Eq. ( ^Tl) are independent of in leading order. 



3.3 Fitting procedure 

In principle, five sets of gi structure function data are available for a global fit: the proton 
data from SLAG- Yale 0, EMC § and SMC fTg, the neutron data from SLAC-E142 §, 
and the deuteron data from SMC p. Our analysis, however, is based on a leading- twist 
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Figure 2: Fit to the gf structure function with set A gluon 
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Figure 3: Prediction for the gf structure function with the set A gluon 
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approximation, and it is therefore necessary to impose a minimum cut. Based on 
similar fits to unpolarized structure function data, see for example |T^, and on recent 
theoretical estimates p^, we believe that > 4 GeV^ is sufficient. This means that 
the bulk of the SLAG data (£'beam = 30 GeV) falls outside our fit region. Only the EMG 
and SMG proton data (-Ebeam ~ 200 GeV) cover a broad range of x for >A GeV^. 

The results of the fit are shown in Figs. ^ and |^. It turns out that the errorbars 
on the deuteron data are so large that there is no statistically significant constraint. 
The fit to Qi, Fig. |], is equally good for the three different types of gluon distribution 
described in Section |0| (see Table ^, and so only the curves corresponding to set A 
are shown. Each data point is fitted at its appropriate x and value, but since data 
from different beam energies are shown on the same figure, the curves correspond to 
fixed = 4^ 10, 30 GeV^ For the SMG and EMG data, these correspond to the small-, 
medium- and large-x data respectively. For the SLAG data, only the largest x data points 
have ~ 4 GeV^. The lower x data points are not used in the fit but are displayed for 
completeness. Figure § shows our predictions for the neutron structure function, together 
with the SLAG-E142 data g. Only the largest x data point (with = 5,2 GeV^) can 
be compared with the theoretical curves. At lower x, we might expect that higher-twist 
contributions are important, although it is interesting that the = A GeV^ curve still 
gives a fair description of the data. 

The values of the fitted parameters and their associated errors are listed in Table |^. 
Not surprisingly, the quark distributions are almost independent of the choice of gluon 
distribution and the quality of the fit is similar for sets A, B and G. Note that the Qu (= dd) 
parameter is reasonably constrained to have a value ~ 0.45, close to the corresponding 
unpolarized values, see 





Gluon A 


Gluon B 


Gluon G 


7g 


0.0 


18.0 


-3.5 


xVdof 


18.8/30 


18.9/30 


20.8/30 




0.46±0.15 


0.45ib0.12 


0.43±0.12 


lu 


18.36±14.49 


20.22±12.2 


17.05±10.71 


be 


7.44±3.52 


11.07±4.19 


8.00±1.59 



Table 4: The fitted parameters and x^/dof values corresponding to the three types of 
gluon distribution. Note that in each case a„ = and 7^ = 7^. 



14 



0.03 
0.02 
0.01 


-0.01 
-0.02 



Figure 4: Prediction for the structure function with the set A gluon 



Figures show the gluon, quark-singlet, m- valence and (i- valence distributions at Ql 
obtained from the fits. For comparison, we show also the unpolarized distributions of 
Ref. [ig. 

A final comment concerns the shape of the gluon distribution at small x. For all of 
our sets, the gluonic contribution to gi is smooth and negative at small x. If one uses a 
different gluon coefficient function, for example either of those in Table |^, then the same 
smooth single-signed behaviour can only be reproduced by using an oscillating gluon, see 
Fig. |l|. Conversely, a smooth gluon combined with either of these alternative coefficient 
functions would give a gluonic contribution to gi which changes sign at small x. 




4 Predictions for /\q 

In our model, a non-zero antiquark distribution Ag is generated at higher by the 
evolution of the quark singlet, Eq. (0). The resulting Ag is sensitive to the initial shape of 
the gluon. A measurement of the polarized sea-quark distribution for > A GeV^ (recall 
that we generate an SU(3) symmetric sea) would therefore probe the gluon distribution. 
There are several ways that this could be done. Perhaps the most promising in the short 
term is to identify specific flavours of hadrons in the final-state. The extension of Eq. (^) 
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Figure 5: The polarized gluon (set A), quark singlet, Uy and dy distributions at Q 
4 GeV^ obtained from the fit to the deep inelastic scattering data. 
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Figure 7: As for Fig. ||, but with the set C gluon distribution. 
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Figure 8: The xAg distributions at = 4, 10, 30 GeV^ corresponding to the (a) set A, 
(b) set B, and (c) set C gluon distributions. 



to hadron inclusive structure functions for Ip — » IH + X is 



9i 



(45) 
(46) 



where D'^^ (z) are fragmentation functions. As explained, for example, in Ref. the 
identification of vr^, K^, ... in the final state allows the sea-quark polarization to be 
probed directly. Our predictions for the polarized sea-quark distribution are shown in 
Fig. |[ Independent of the choice of gluon, all the distributions have the common feature 
of being negative for small x < (9(0.01). At large x, the fact that AGc is less than zero 
(see Fig. |^) is reflected in the corresponding Aqc distribution. 

A preliminary study of the final-state hadron charge asymmetry in polarized deep 
inelastic scattering by SMC has given a first indication of the size and shape of the Aq 
distribution. Although the experimental errors are still rather large, all the distributions 
of Fig. P are consistent with these data. 



5 Summary and conclusions 

In this paper we have described an attempt to perform a global fit to polarized deep 
inelastic structure function data on proton targets in order to extract a consistent set 
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of polarized parton distributions. We have adopted the point of view that there is no 
significant polarization of the sea quarks at low Q^, and therefore that the gluon carries a 
substantial fraction of the proton's spin. We have constrained the net spin carried by the 
various parton flavours using deep inelastic sum-rule and hyperon decay measurements. 
Recognizing that the quality of the data is not yet high enough to fully constrain the 
shapes of the distributions, we have imposed reasonable theoretical constraints at large 
and small x. We have performed fits with three qualitatively different gluon distributions, 
chosen to reflect the extent to which fast-moving gluons 'remember' the spin of the parent 
proton. 

Our distributions have several applications: 

(i) they provide a benchmark for future deep inelastic structure function measurements 
- in particular, measurements using neutron and deuteron targets, 

(ii) they can be compared with theoretical, non-perturbative model calculations of the 
spin structure of the proton, 

(iii) they can be used to predict cross sections for the production of various flavours 
of hadrons in the final state - vr^ and K"^ for probing the quark sea, and the 
production of heavy quarks (in particular J /if) and open charm) for probing the 
gluon distribution, and 

(iv) they can be used as input for the calculation of hard scattering cross sections in 
polarized hadron-hadron collisions. 

In conclusion, polarized lepton-hadron scattering experiments are beginning to pro- 
vide precision information on how the spin of the proton is shared among its parton 
constituents. This will be invaluable in testing the predictions of Quantum Chromody- 
namics and in designing future experiments. Our distributions should be a useful tool in 
both these areas. Q 
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